Evolution of a collapsing and exploding Bose-Einstein condensate in different trap 

symmetries 
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Based on the time-dependent Gross-Pitaevskii equation we study the evolution of a collapsing and 
exploding Bose-Einstein condensate in different trap symmetries to see the effect of confinement on 
collapse and subsequent explosion, which can be verified in future experiments. We make prediction 
for the evolution of the shape of the condensate and the number of atoms in it for different trap 
symmetries (cigar to pancake) as well as in the presence of an optical lattice potential. We also 
make prediction for the jet formation in different cases when the collapse is suddenly terminated by 
changing the scattering length to zero via a Feshbach resonance. 
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I. INTRODUCTION 

Since the detection and study of Bose-Einstein conden- 
sates (BECs) of ''Li atoms with attractive interaction 0, 
such condensates have been used in the study of solitons 
and collapse Q- In general an attractive condensate 
with number of atoms N larger than a critical value A^cr 
is not dynamically stable 0. However, if such a strongly 
attractive condensate is "prepared" or somehow made 
to exist it experiences a dramatic collapse and explodes 
emitting atoms. The first demonstration of such a col- 
lapse was made with a ''Li condensate by slowly increas- 
ing the number of atoms in it from an external source, 
while the BEG showed a sequence of collapse with the 
number of atoms N oscillating around iVcr- Such a col- 
lapse is driven by a stocastic process. 

A dynamical study of a much stronger and violent col- 
lapse has been performed by Donley et al. Q on an 
attractive ^^Rb BEG Q in an axially symmetric trap, 
where they manipulated the inter-atomic interaction by 
changing the external magnetic field exploiting a nearby 
Feshbach resonance |5J. In the vicinity of a Feshbach 
resonance the atomic scattering length a can be varied 
over a huge range by adjusting the external magnetic 
field. Gonsequently, they changed the sign of the scatter- 
ing length, thus transforming a repulsive condensate of 
^^Rb atoms into an attractive one which naturally evolves 
into a collapsing and exploding condensate. Donley et al. 
provided a quantitative estimate of the explosion of this 
BEG by measuring different properties of the exploding 
condensate. 

It has been realized that many features of the exper- 
iment by Donley et al. M on the collapsing condensate 
can be described H 0, ES IH 111 111 111 111 III 
m Eiby the mean-field Gross-Pitaevskii (GP) equa- 
tion [13. To account for the loss of atoms from the 
strongly attractive collapsing condensate an absorptive 
nonlinear three-body recombination term is included in 
the GP equation 0. However, we are fully aware that 
there are features of this experiment which are expected 



to be beyond mean-field description. Among these are 
the distribution of number and energy of emitted high- 
energy ('^ I0~^ Kelvin) uncondensed burst atoms re- 
ported in the ex peri ment. Although there have been 
some attempts [3, 111 Ell to describe the burst atoms us- 
ing the mean-field GP equation, now there seems to be a 
consensus that they cannot be described adequately and 
satisfactorily using a mean- field approach jl^, El Il5| . 
Also, the GP equation does not successfully predict the 
"time to collapse" (or the time lag to start the collapse af- 
ter changing the sign of the scattering length) in all cases 
investigated in the experiment, as has been pointed out 
inRefs. El El. 

The GP equation is supposed to deal with the zero- 
or very low-energy condensed phase of atoms and has 
been used to predict the time to collapse, evolution of 
the collapsing condensate including the very low-energy 
(~ nano Kelvin) jet formation H when the collapse is 
suddenly stopped before completion by jumping the scat- 
tering length to Oquench ~ ^ (noninteracting atoms) 
or positive (repulsive atoms) values. The jet atoms are 
slowly formed in the radial direction when the collapse 
is stopped in this fashion. In the experiment usually 
'^quench = 0. It is emphasized that unlike the emit- 
ted uncondensed "hotter" missing and burst atoms re- 
ported in the experiment I3| the jet atoms form a part 
of the surviving "colder" condensate and hence should 
be describable by the mean-field GP equation. Saito et 
al. 0, Bao et al. El ^ind this author E3 presented a 
mean- field description of jet formation and Galzetta et al. 
|l4j treated jet formation exclusively as a quantum effect. 
More recently, the present author has used a set of cou- 
pled mean-field-hydrodynamic equations El to describe 
the essentials of the collapse dynamics of a mixture of a 
boson and fermion condensates |20l |. 

In this paper we extend the study of the evolution of 
the collapsing and exploding condensate in different sym- 
metries to see the effect of confinement on collapse and 
subsequent explosion. Future experiments may verify 
these predictions and thus provide a more stringent test 
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for the mean-field GP equation. The experiment of Don- 
ley et al. was performed for an axially-symmetric cigar- 
shaped BEG. In the present analysis we extend our study 
to a spherical as well as an axially-symmetric pancake- 
shaped BEG. 

Lately, the periodic optical-lattice potential has played 
an essential role in many theoretical and experimental 
studies of Bose-Einstein condensation, e. g., in the study 
of Josephson oscillation l^^l and its disruption , inter- 
ference of matter- wave |23j, BEG dynamics on periodic 
trap etc. The optical-lattice confinement creates a 
BEG in an entirely different shape and trapping condi- 
tion form a conventional harmonic oscillator trapping. 
Gonsequently, one could have a collapse of a different na- 
ture in the presence of an optical-lattice potential. We 
shall see in our study that under certain conditions of 
trap symmetry, in addition to the usual global collapse 
to the centre, in the presence of the optical-lattice poten- 
tial one could have independent local collapse of pieces 
of the condensate to local centres. In view of this we 
study the dynamics of a collapsing and exploding BEG 
of different symmetries prepared on a periodic optical- 
lattice potential. We study the evolution of the shape 
and size of the condensate as well as the jet formation 
upon stopping the collapse by making the BEG repulsive 
or noninteracting. 

In Sec. II we present our mean-field model. In Sec. Ill 
we present our results that we compare with the exper- 
iment and other numerical studies. In Sec. Ill we also 
present a physical discussion of our findings and some 
concluding remarks are given in Sec. IV. 



the axial direction z of the harmonic trap. The cigar- 
shaped condensate corresponds to v < 1 and pancake- 
shaped condensate corresponds to v > 1. The periodic 
optical-lattice potential in the axial z direction created 
by a standing- wave laser field of wave length A is given by 
V^p = KEiiCos^{kLz) with Er — 1i^k\/ {2m),kL = 2tt/\ 
and K the strength. The normalization condition of the 
wave function is / (ir|\I'(r; r)p = 1. Here we simulate 
the atom loss via the most important quintic three-body 
term if 3 0, 0, H, • The contribution of the cubic two- 
body loss term K2 [l^ is expected to be negligible 0, Q 
compared to the three-body term in the present problem 
of the collapsed condensate with large density and will 
not be considered here. 

In the absence of angular momentum the wave function 
has the form '^{r\ r) = il){p, z; t). Now transforming to di- 
mensionless variables defined by a; = V2p/l, y — ^/2z/l, 
t — TOJ, I = y^h/ (mw), and 



0(x,2/;O ^ ^i^lM = J-^^{p,z;r), (2.2) 
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II. NONLINEAR GROSS-PITAEVSKII 
EQUATION 

The time-dependent Bose-Einstein condensate wave 
function 5'(r;T) at position r and time r allowing for 
atomic loss may be described by the following mean-field 
nonlinear GP equation 



- + + gAf|*(r;T)|2 - - 



X {K2N\^{r;T)\^ +K3N^\'f{r;T)\'' 



^'(r;T) = 0. 

(2.1) 



Here m is the mass and N the number of atoms in 
the condensate, g = Anh^a/m the strength of inter- 
atomic interaction, with a the atomic scattering length. 
The terms K2 and K3 denote two-body dipolar and 
three-body recombination loss-rate coefficients, respec- 
tively and include the Bose statistical factors 1/2! and 
1/3! needed to describe the condensate. The trap po- 
tential with cylindrical symmetry may be written as 
V{r) = ima;^(p^ -I- z/^z^) -I- Vbp where uj is the angu- 
lar frequency in the radial direction r and vlu that in 



where n = Na/l, Xq = V2\/l and C = iKs/^a'^l^u;). The 
normalization condition of the wave function becomes 



= 2tt 



dx 



dy\ip{x,y;t)\^x-\ (2.4) 



For C — ~ 0, Wnorm = 1, however, in the presence 
of loss K3, > 0, A/norm < 1- The number of remaining 
atoms N in the condensate is given by A'^ = A^oA/norm, 
where Nq is the initial number of atoms. 

In this study the term K3 or — 4X3(0^1*11!) will be 
used for a description of atom loss in the case of attrac- 
tive interaction. The choice of has a huge effect on 
some experimental observables and the fact that it is ex- 
perimentally not precisely determined is a problem for 
existing theory on the experiment. As in our previous 
study 'ITI we employ C = 2 and K3 ^ throughout this 
study. It was found QJ^ that this value of C(= 2) repro- 
duced the time evolution of the condensate in the exper- 
iment of Donley et al. ^ satisfactorily for a wide range 
of variation of initial number of atoms and scattering 
lengths |1. The present value C = 2 with K3 = 
leads to ifg ~ 8 x IQ-^^ cmVs at a = -340ao and 
~ 6 X 10"'^^ cm^/s at a = — 30ao, where uq is the 
Bohr radius. The experimental value of loss rate is 
A'3 ~ 7 X 10~^^ cm^/s at a = — 340ao which is very close 
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to the present choice. Of the theoretical studies, the 
values used by Santos et al. |l3 [K^ ~ 7 x 10~^^ cm^/s 
at a = -340ao), Savage et al. [ij (i^s ~ 19 x lO^^T 
cm^/s at a = -30ao), Bao et al. [isj (i^a ^ 6.75 x IQ-^"^ 
cm^/s at a = — 30ao) and the present author are con- 
sistent with each other and describes well the decay of 
the collapsing condensate. 



III. NUMERICAL RESULT 
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FIG. 2: (Color online) The time to collapse iggllapse ^® 
I '^collapse I for Ojjjj^jg^j = 0, A^o = 6000. Solid circle with 
error bar: experiment Hwith v = 0.39 and k = 0; Open cir- 
cle: mean-field model of Q with v = 0.39 and /t = 0; full line: 
present result with v = 0.39 and ti = Q; dashed line: present 
result with u = 1 and k = 0; dashed-dotted line: present re- 
sult with I' = 0.1 and «: = 4; and dashed-doubled-dotted line: 
present result with = 0.1 and k = 0. 



FIG. 1: (Color online) Number of remaining atoms N{t) in the 
condensate of 16000 ^^"'Rb atoms after ramping the scattering 
length from ajj^ = 7ao to (a) a^gj = — 30ao and (b) a^^j = 
— 6.7ao as a function of evolution time in milliseconds. The 
unpublished and unanalyzed experimental points of Donley et 
al. |3| for a^Qj = — 6.7ao are taken from Bao et al. d. The 
curves are labeled by their respective optical lattice strength 
K and axial trap parameter v. 

We solve the GP equation 1)2.3(1 numerically using a 
time-iteration method based on the Crank-Nicholson dis- 
cretization scheme elaborated in |2^. We discretize the 
GP equation using time step A = 0.001 and space step 
0.1 for both X and y spanning x from to 15 and y from 
—30 to 30. This domain of space was sufficient to encom- 
pass the whole condensate wave function in this study. 



First, the numerical simulation is performed with the 
actual parameters of the experiment by Donley et al. , 
e. g., the initial number of atoms, scattering lengths, 
etc. Throughout this investigation we take the harmonic 
oscillator length I = 2607 nm and one unit of time t = 
0.009 095 s i3 consistent with the experiment of Donley 
et al. Q. When we include an optical-lattice potential, 
the optical-lattice strength n is taken to be 4, and the 
reduced wave length Aq is taken to be 1 throughout this 
study. These optical-lattice parameters are consistent 
with the experiment by Gataliotti et al. 0, |23| ■ The 
numerical simulation using Eq. 1(2.3(1 with a nonzero C(= 
2) immediately yields the remaining number of atoms in 
the condensate after the jump in scattering length. 
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A. Evolution of the Number of Atoms in the 
Condensate 

In the experiment the initial scattering length a-^j^(> 0) 
of a repulsive condensate is suddenly jumped to 0'(.Qii< 0) 
to start the collapse. The remaining number N{t) of 
atoms vs. time for an initial number of atoms Nq = 16000 
and an initial scattering length ajj^ = Tag are shown in 
Figs. 1 (a) and (b) for final scattering lengths after col- 
lapse a^^pj = — SOflQ and —6. Tag, respectively. In both 
cases the experimental data for k = and v = 0.39 
(cigar-shaped condensate) are in agreement with the the- 
oretical simulation without any adjustable parameter. 
For Oj^qJ = — 6.7ao, the unpublished experimental data 
of 3| as shown in Fig. 1 (b) are as quoted in Bao et 
al (ISj]. These data are not fully analyzed and for large 
time are expected to be bigger than the actual number 
of atoms. This is due to the difficulty in separating the 
remnant condensate from the oscillating atom cloud sur- 
rounding it 3]. In addition, in Figs. 1 we plot the results 
for K = 4 and v = 0.39 (cigar-shaped condensate with 
optical-lattice potential) ; n — A and v — 1 (spherical con- 
densate with optical-lattice potential); k = and v — 1 
(spherical condensate); and k = and v = b (pancake- 
shaped condensate). 

As the repulsive condensate is quickly turned attrac- 
tive at <: = 0, via a Feshbach resonance, the condensate 
starts to collapse and once the central density increases 
sufficiently it loses a significant portion of atoms in an 
explosive fashion via three-body recombination to form 
a remnant condensate in about 15 ms as can be seen in 
Figs. 1. After explosion the number of atoms in the 
remnant continues to be much larger than the critical 
number of atoms A'cr and it keeps on losing atoms at a 
much slower rate without undergoing violent explosion. 
However, in some cases the remnant undergoes a smaller 
secondary explosion while it loses a reasonable fraction 
of atoms in a small interval of time. This happens when 
the number of atoms in the remnant is much larger than 
Nqy so as to initiate a secondary collapse and explo- 
sion. Prominent secondary explosions in the presence 
of optical-lattice potential are found in different cases in 
Figs. 1 for 40 > t > 30. 

B. Time to Collapse 



for a cigar-shaped symmetry. These features of time to 
collapse are illustrated in Fig. 2 where we plot ^eoUapse 
vs. kcollapsel/^o collapse of a condensate of 6000 

atoms originally in a noninteracting state with scattering 
length ajjiii^jai — 0. Then suddenly its scattering length 
is changed to a negative (attractive) value ^collapse ^'^^ 
its i(2ollapse obtained. Donley et al. experimentally 
measured ^collapse ^^^^ ''^^^ ^^'^ ^ ~ ^'^^ ^^"^ k = 
and here we provide the same for other values of trap 
symmetry v and also in the presence of a optical lat- 
tice potential with k — A. It should be recalled that the 
prediction of the GP equation by this author and others 
[3 [T3 . IT^ does not very well describe the experimental 
results of Donley et al. for the time to collapse. The 
inclusion of the optical-lattice potential has reduced the 
time to collapse in a cigar shaped condensate {v — 0.1). 

The above features of time to collapse could be un- 
derstood on a physical ground. In a cigar-shaped con- 
densate the average distance among the atoms is larger 
than that in a pancake-shaped condensate of same vol- 
ume. Hence, due to atomic attraction a cigar-shaped 
condensate has to contract during a larger interval of 
time than a pancake-shaped condensate before the cen- 
tral density increases sufficiently to start an explosion. 
This justifies a larger time to collapse for a cigar-shaped 
condensate. In the presence of an optical-lattice potential 
for cigar-shaped symmetry the optical-lattice divides the 
condensate in a large number of pieces. What predomi- 
nates in the collpase of such a condensate is the collapse 
of an individual piece to a local center rather than to 
the global center of the condensate via tunneling. This 
is a quicker process than the collapse of the whole con- 
densate to the global center. This is why the time to 
collpase is shorter for a cigar-shaped condensate in an 
optical-lattice trap than a cigar-shaped condensate in a 
harmonic trap alone. In a pancake-shaped symmetry the 
number of optical-lattice sites inside the condensate is 
small. In this case a separation of the condensate in a 
smaller number of pieces does not aid in the collapse, as 
the different slices of the condensate has to collapse es- 
sentially towards the center of the condensate before the 
explosion starts. Hence the optical-lattice potential has 
almost no effect on the time to collapse in the pancake- 
shaped or spherical symmetry. 



Another important aspect of collapse is the "time to 
collapse" or the time to initiate the collapse and explo- 
sion ^collapse ^ft^'" the repulsive condensate is suddenly 
made attractive ai t — Q. Collapse is characterized by 
a sudden rapid emission of atoms from the condensate. 
From Figs. 1 we find that the time to collapse is the 
shortest for a pancake-shaped symmetry {v > 1) and is 
the longest for a cigar-shaped symmetry [v < 1). The in- 
clusion of an optical-lattice potential has no effect on the 
time to collapse for a spherical or pancake-shaped sym- 
metry. However, its inclusion reduces the time to collapse 



Another aspect of Figs. 1 worth mentioning is that the 
number of atoms in the remnant condensate after the first 
explosion is larger in the presence of an optical-lattice 
potential. Due to optical-lattice barriers one essentially 
has local collapse of different pieces of the condensate in 
this case as opposed to a global collapse to the center 
of the condensate in the case of a harmonic trap alone. 
Consequently, the collapse is more violent with greater 
loss of atoms in the absence of an optical-lattice trap. 
This is why the remnant number after the first explosion 
is larger in the presence of an optical trap. 
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(a) K = 0, v = 0.39 




(b) K = 4, v = 0.39 




FIG. 3: A view of the evolution of the residual condensate 
wave function \tp{p, z)\ in arbitrary units for initial scattering 
length ajjj — 7ao, final scattering length a^^y — — 30ao, initial 
number of atoms A^o ~ 16000 at times t — (i) 4 ms, (ii) 6 
ms, (iii) 8 ms and (iv) 12 ms for (a) k = 0, = 0.39 and (b) 
K = 4, = 0.39. 



C. Evolution of the Shape of the Condensate 

Next we consider the evolution of the shape of the 
residual condensate. In Fig. 3 (a) we show the profile 
of the wave function ip^p, z) at different times during ex- 
plosion for iVo = 16000, K = 0, = 0.39, a^^ — Tuq, and 
Aj^qJ = — 30ao. This is the case of a cigar-shaped con- 
densate used in the experiment j3| . During explosion the 
condensate wave function develops a three-peak structure 
noted before in P| . In Fig. 3 (b) we illustrate the profile 
of the wave function z)| at different times during 
explosion of the condensate formed in an optical-lattice 
potential with k = 4 in addition to the axial harmonic 
trap: other parameters remaining the same as in Fig. 3 
(a). The condensate now develops a distinct multi-peak 
structure along the optical lattice in place of the three- 




FIG. 4: Same as in Figs. 3 for ft = 0, = 5. 

peak structure in the absence of the optical-lattice poten- 
tial. However, the number of peaks in the wave function 
is less by a factor of two to three than the number of pits 
of the optical-lattice potential. The number of distinct 
peaks in the wave function in this case is five as can be 
seen in Fig. 3 (b). 

The above distinct peaks in the wave function in the 
presence of the optical-lattice potential may have inter- 
esting application in the generation of radially bound and 
axially free bright solitons. The wave function of Fig. 3 
(b) is axially bound. However, if the axial trap and the 
optical-lattice potential are removed, or better an expul- 
sive potential is applied in the axial direction, the wave 
function will expand axially. The side peaks of the wave 
function can evolve into separate solitons and come out 
in the axial direction which can be used as bright solitons 
in other experiments. 

The scenario of the evolution of the condensate is 
entirely different for pancake-shaped condensate with 
v > 1. In that case the condensate is squezeed in the 
axial direction and a single peak, rather than multiple 
peaks, is formed in the condensate wave function. This 
is illustrated in Fig. 4 where we plot the condensate 
wave function for k = and v = 5, the other param- 
eters of simulation being the same as in Figs. 2. The 
use of optical-lattice potential in this case also does not 
lead to prominent peaks in the wave function in the axial 
direction. 



D. Jet Formation 

Another interesting feature of the experiment of Don- 
ley et al [3 is the formation of jet. As the collapse was 
suddenly terminated after an evolution time tev by jump- 
ing the scattering length from a^^-^ to o-q^Quch. — ^' 
jet atoms were slowly formed in the radial direction. In 
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FIG. 5: (Color online) A view of the evolution of radial jet 
at times t — 0, 2 ms, 4 ms and 5.2 ms on a mat of size 
16 fim X 16 fim from a contour plot of |?/'(p, z)] for initial 
scattering length a^j^ — 7ao, final scattering length a^^^ ~ 
— 30ao, initial number of atoms A^o — 16000, (a) without an 
optical-lattice potential (k — 0) and (b) with an optical-lattice 
potential (k = 4). In both cases the jet formation was started 
by jumping the scattering length to iquench ~ ^ after a time 
iev = 4 ms of the beginning of collapse. 



the strongly collapsing condensate, local radial spikes are 
formed during particle loss as can be seen from a plot of 
the numerically calculated wave function and in ex- 
periment 3] . During particle loss the top of the spikes 
are torn and ejected out and new spikes are formed un- 
til the explosion and particle loss are over. There is a 
balance between central atomic attractive force and the 
outward kinetic pressure. If the attractive force is now 
suddenly removed by stopping the collapse by applying 
"^quench ^ ^^'^ highly collapsed condensate expands 
due to kinetic pressure, becomes larger and the recom- 
bination of atoms is greatly reduced. Consequently, the 
spikes expand and develop into a prominent jet 0. 

Now we consider the jet formation as in the experiment 
of Donley et al. Q at different times t of the collapsing 



FIG. 6: (Color online) Same as in Figs. 5 for a^^-^ = — 250ao, 
and tev ~ 2 ms. 




FIG. 7: (Color online) Same as in Fig. 5 (a) for v = 4 and 
t = 0, 2 ms, 4 ms and 6 ms. 
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condensate after jumping the scattering length suddenly 
from a^^Y = — 30ao to aq^g^cj-^ = during explosion at 
time tev from the beginning of collapse and explosion. 
The initial scattering length aj^ = 7ao and number of 
atoms A^o = 16000. In Figs. 5 (a) and (b) we show the 
contour plot of the condensate for iev = 2 ms, without 
(k = 0) and with (k = 4) an optical-lattice potential, 
respectively, at different times t = 0, 2 ms, 4 ms, and 5.2 
ms after jumping the scattering length to iquench ~ ^' 
A prominent radial jet is formed slowly at time t = 4 — 6 
ms after stopping the collapse at tev = 4 ms. The jet 
is much less pronounced for tev — 2, 8 ms and 10 ms 
(not shown in figure) compared to the jet in Figs. 5. 
There is a fundamental difference between the jets in 
Figs. 5 (a) and (b) in the absence and presence of optical- 
lattice potential. In Fig. 5 (a) the absence of the optical 
potential the jet is narrow, whereas in Fig. 5 (b) it is 
wide and spread over a number of optical lattice sites. 

In addition, we studied jet formation for different val- 
ues of cl^qY i'^ place of a^j^j = — 30ao and find that 
the general scenario remains similar. For example, for 
(IqqI = — 250ao, the collapse and subsequent explosion 
starts at a small time. So for a good formation of jet 
a smaller value of iev is to be preferred. In Figs. 6 we 
show the jet formation for a^^Y = — 250ao and tev — 2 
ms. In this case the shape of the jet is different from that 
in Figs. 5. However, as in Figs. 5, the jet gets broadened 
in the presence of the optical-lattice potential. 

Next we study the effect of the axial trap symmetry 
on jet formation. In Figs. 5 and 6 the harmonic trap 
has cigar symmetry {v — 0.39 < 1). If it is changed to 
pancake symmetry (v > 1), the condensate and the jet 
gets compressed in the axial direction. Consequently, the 
radial jet is not very pronounced and can not be clearly 
distinguished from the condensate. This is illustrated in 
Fig. 7 for pancake-shaped trap with v = Ain the absence 
of an optical-lattice potential. In this case at t = 6 ms 
the condensate is more extended in the radial direction 
compared to the condensate at t = due to the formation 
of jet. However, due to the overall compression of the 
condensate in the axial direction the jet can not be easily 
separated from the condensate. In contrast in Figs. 5 the 
jet is easily separated from the central condensate. So 
pancake symmetry is not ideal for the study of a jet. The 
situation does not change in the presence of an optical- 
lattice potential superposed on a pancake-shaped trap. 

IV. CONCLUSION 

In conclusion, we have employed a numerical simula- 
tion based on the accurate solution "2^ of the mean-field 



Gross-Pitaevskii equation with a cylindrical trap to study 
the evolution of a collapsing and exploding condensate 
as in the experiment of Donley et al. In the GP 

equation we include a quintic three-body nonlinear re- 
combination loss term 6] that accounts for the decay of 
the strongly attractive condensate. We also extend our 
investigation to different trap symmetries and including 
an optical-lattice potential in the axial direction. In ad- 
dition to studying the evolution of the size and the shape 
of the condensate, we also study the jet formation as ob- 
served experimentally. Without any adjustable parame- 
ter the result of the present and previous simulations of 
this author are in good agreement with some aspects of 
the experiment by Donley et al. 0- 

It is interesting to emphasize that the GP equation 
does describe some but not all aspects of the collapse ex- 
periment by Donley et al. and its predictions for the 
"time to collapse" do not agree well with experiment 
[T^ IT^ . The failure to explain time to collapse is dra- 
matic as intuitively one should expect the mean-field 
model to be a faithful model for time to collapse involv- 
ing the dynamics of the coldest atoms in the condensate. 
However, there are aspects of experiments which cannot 
be described by mean- field models 0, IT^ H^IT^. e. 
g., the dynamics of missing and burst atoms [Sj. Fur- 
thermore all numerical studies of this experiment suffer 
from limited knowledge of the three-body loss rate 
and even though many experimental features can be de- 
scribed by a suitable choice of K^, no value of K3 yields 
simultaneous agreement between predictions of the GP 
equation and all observable quantities of the experiment. 
In this situation it would be of great theoretical and ex- 
perimental interests to see if a repeated experiment with 
different trap parameters (and also including an optical- 
lattice potential) would help to understand the underly- 
ing physics and could make use of numerical studies such 
as those described in the present paper. 
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